MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations

Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain statistically exhaustive models of the long-timescale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large lengthscales, such as in cellular signal transduction, particle-based reaction-diffusion (RD) simulations are more suitable but come at the price of sacrificing molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large timeand lengthscales, while still conserving the characteristic features of the interactions observed at atomic detail. While the insertion of a MSM rate matrix into an RD framework seems formally straightforward, fundamental questions about the coupling are open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make a first step towards MSM/RD by laying out a general theory of coupling and proposing a first implementation that includes both conformational changes (A B) and association/dissociation of proteins with small ligands or substrates (A + B C). Applications on a toy model and CO diffusion and binding to myoglobin are reported. †Equal contribution 1Freie Universität Berlin, Department of Mathematics and Computer Science, Arnimallee 6, 14195 Berlin, Germany 2CIC nanoGune, San Sebastián 20018, Spain a)Corresponding author. Electronic mail: [email protected].